Quantum phase transitions in the quasi-periodic kicked rotor 
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We present a microscopic theory of transport in quasi-periodically driven environments ('kicked 
rotors'), as realized in recent atom optic experiments. We find that the behavior of these systems 
depends sensitively on the value of Planck's constant h: for irrational values of /i/(47r) they fall 
into the universality class of disordered electronic systems and we derive the microscopic theory 
of the ensuing localization phenomena. In contrast, for rational values the rotor-Anderson insu- 
lator acquires an infinite (static) conductivity and turns into a 'super-metal'. Signatures of the 
corresponding metal/super-metal transition are discussed. 
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The quantum kicked rotor (QKR) is one of the most 
prominent model systems of quantum nonlinear dynam- 
ics (quantum chaos) 1 . It comprises a one-dimensional 
quantum particle subject to periodic boundary condi- 
tions and sequential driving in time. In spite of the nom- 
inal simplicity of this system — the QKR Hamiltonian 
is defined in terms of only two dimensionless parameters, 
kicking strength, K, and time in relation to (effective) 
Planck's constant, ft. [2J — it displays striking parallels 
to disordered electronic systems [31 H] , and the ensuing 
connections to condensed matter physics have been a sub- 
ject of fundamental research for more than three decades 
(cf. Ref. |4| and references therein). The general inter- 
est in the QKR took a leap forward when the system was 
experimentally realized in a cold atom setting [S| , and lo- 
calization phenomena otherwise predicted for disordered 
multi-channel quantum wires were observed. This has 
been the first observation of strong Anderson localiza- 
tion in a quasi-one dimensional setting (although the full 
microscopic correspondence to disordered quantum wires 
was established only recently f?'.) The next experimental 
breakthrough occurred in 2008 when an effectively higher 
dimensional rotor, the so-called quasi-periodic kicked ro- 
tor, was realized in a gas of cold cesium atoms, and an 
Anderson transition was observed [6 . By now, even sig- 
natures of the critical states emerging at the transition 
point have been seen 'Tj , and it stands to reason that the 
quasi-periodic QKR makes for an almost ideal environ- 
ment to study Anderson type critical phenomena. 

At the same time, the correspondence between the 
quasi-periodic QKR and {d > l)-dimensional disordered 
systems is not as well understood as in the d = 1 case, and 
this is a gap which we aim to close in this Letter. In the 
quasi-periodic QKR, d-dimensional behavior is simulated 
by modulated driving at d different frequencies. Below, 
we will map the low energy physics of this system onto an 
effective field theory equivalent to the non-linear cr-model 
of disordered metallic systems [Sj. For generic values of 
the system parameters, this construction establishes the 
connections to disordered systems, and it explains the 
observation of Anderson type criticality. However, the 



rotor is not a genuine metal, and these differences show 
in anomalies at certain (non-generic) parameter values. 
Specifically, the one-dimensional QKR displays so-called 
quantum resonances (cf. Ref. [4 and references therein) 
at rational values h = Anp/q, p,q (zN coprime. At res- 
onance, the system behaves like a quasi-one dimensional 
ring (in the space of rotor-angular momenta) of radius q, 
and no Anderson localization is observed. 
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FIG. 1. The d = 2 QKR at resonance. Left: g-periodic an- 
gular momentum space. Right: compactification to cylinder 
of circumference q, threaded by flux <j). Grey shaded areas 
represent a 2d lattice of angular momentum sites. Transport 
between sites A'^ and N' takes place through phase coherent 
propagation of advanced (solid) and retarded (dashed) quan- 
tum amplitudes subject to different flux values (l)±. Bottom, 
the three distinct regimes in virtual directions, a) unbound 
diffusion, b) localization at ^ < g, and c) localization at ^ > g. 



Below we will show that in dimensions d > 1 the sys- 
tem displays even richer behavior. At ft = Anp/q the 
system compactifies, e.g., for d = 2 a nominally two di- 
mensional QKR maps onto the surface of an infinitely 
long cylinder of finite circumference q (cf . Fig. [l] top) . 
As far as bulk localization phenomena are concerned, this 
entails a dimensional reduction to quasi d~l dimensions 
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(localization along the cylinder axis). However, measur- 
able response functions dominantly couple to its compact 
dimension, and this has a number of interesting conse- 
quences: below we consider a (measurable) observable 
which, in the light of the above analogies to metallic sys- 
tem, plays a role analogous to an optical conductivity. 
We find that in the quasi {d — l)-dimensional localized 
regime, the system actually shows 'super-metallic' be- 
havior (diverging conductivity). In contrast, for quasi 
(d— l)-dimensional metallic regimes (e.g., the quasi three- 
dimensional metallic phase realized by driving at four dif- 
ferent frequencies), the conductivity remains finite. For 
d > 3, an Anderson transition in quasi (d— l)-dimensions 
may be driven by changing the kicking strength, and this 
manifests itself in a metal/super- metal transition in ob- 
servable quantities. A survey of the three different phases 
realizable in the QKR is shown in table [l] 



properties in quasi 
(d — l)-dimensions 


non-generic Planck's 
constant 


generic Planck's 
constant 


localized 


super-metallic, 
infinite conductivity 


insulating, 

zero conductivity 


metallic 


metallic, (d > 4) 
finite conductivity 


metallic, (d > 3) 
finite conductivity 



TABLE I. Phases realizable in the quasi-periodic QKR 



In dimensionless units [2 , the time dependent Hamil- 
tonian of the system is defined as H{t) — ^h'^fi^ + 
KV{t) J2m ^(^ ~ ™)- Here, 9 is the rotor's angular vari- 
able, n ~ —idg the angular momentum operator, and 
V{t) = y(cos0, cos(6'i + LJit), ...,cos{6d-i +ujd-it)), the 
kicking potential. The (d — 1) frequencies, uji, ...,ujd~i, 
are incommensurate to the kicking frequency 27r and 
among themselves, and 9i, ...,9d-i are d — 1 arbitrary 
phases. We assume V to be symmetric with respect to its 
d arguments, and of unit characteristic variation. Later 
we will see that the detailed form of V determines the 
diffusion constant of the system, but is largely inessen- 
tial otherwise. We assume h/Air — p/q where the limit of 
an irrational value may be taken by sending p,q ^ oo. 

As in the experimental applications, we consider the 
evolution of a wave function initially uniform in 9. The 
localization properties of such states are probed by the 
correlation function 



i?(t)^ 2ElHn^=i?^'M|0)Pn^ (1) 
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where U'{m) = exp^i^n^) exp{ij-V{m)) is the Flo- 
quet operator and the overline stands for the average 
over the parameters 9i. The mapping to an effec- 
tively higher dimensional system [9] is achieved by in- 
terpreting 1 9) = \9,9i,. . . ,9d-i) as a d-dimensional co- 
ordinate vector, comprising a 'real' angular coordinate, 
and a generalization of the parameters 9i to 'virtual' 
coordinates. Similarly, we introduce a d-dimensional 



angular momentum state, \N) = \n,ni, . . . ,nd-i) 
where ni = —ide- is conjugate to 9i, with eigen- 
values rii G Z. Defining the operator $(m) = 
exp(—im'^^ujihi), it is then not difficult to verify 
that the 'gauge transformed' Floquet operator U{m) = 
$(r7i + 1)U' {m)^~^ (m) becomes time-independent, viz. 
U = f(N)W{Q), f{N) = e'l"'+»E.".ft>^ ^(e) = 
exp[^V{cos9,cos9i,- ■ ■ ,cos9d-i)], and that E{t) = 

5 Eat l(^|t^*|0)P'^^- We have, thus, traded the time 
dependence of the original problem for an effective ex- 
tension to a multi-dimensional Hilbert space spanned by 
the states \N). 

The effective Floquet operator U possesses two fun- 
damental symmetries: time reversal symmetry |10| T : 
t — )> —t, O — ?> —Q, and invariance under the translation 
n ^ h + q. Exploiting the latter, it is straightforward to 
verify that the variable E{t) affords the representation 

27r/q 


where 'tr' is a trace over all states |A^) whose real co- 
ordinate n e {0, ■ ■ ■ ,q} is restricted to a compact 'unit 
cell', and we have defined the 'Bloch-Floquet' operator 
= fiN)W{Q + (l)), where Q + cf) = {9+(t), 9i,..., 6'<j„i). 
As in the quantum mechanics of periodic structures, the 
summation over a Bloch phase, cf), enables us to com- 
pactify n-space to a ring of circumference q with periodic 
boundary conditions. The shift of the angular variable 
9 ^ 9 + <j) shows that (j> couples to the system as an 
Aharonov-Bohm flux, cf. Fig. [T] top. 

Eqs. ^ and ^ are different (yet equivalent) ways 
of probing the time dependent spreading of angular mo- 
mentum states. Anticipating a competition of classical 
diffusion and quantum localization, we expect three qual- 
itatively distinct cases (cf. Fig. [l] bottom): if the lo- 
calization length, ^ is infinitely large, unbound diffusive 
spreading ~ Dt characterized by a diffusion coeffi- 
cient, D, leads to a linear increase E{t) ~ Dt {metallic 
regime). In contrast, for ^ < q we expect saturation, 
E{t ^ ^'^/D) ^ const, {localized regime). (A finite size 
correction ~ exponentially small in e~'^^^ will be ig- 
nored.) Finally, in cases where > q, the system be- 
haves similar to a finite quantum system of character- 
istic quasilevel spacing ~ l/{q^''-~^). For large times, 
t ^ q£,^-'''~'^\ individual states of this system can be re- 
solved and a formal decomposition of U in quasi-energy 
states shows that E{t) ~ {super-metallic regime). 

To quantitatively describe these regimes, we define 
the resolvent operators GJ(w±) = (1 - (e''^±f70)±i)-\ 
where uj± = (jJq zt ^{uj -\- iO). The Fourier transform 
E{u}) = dte^'^*E{t), then assumes the form E{uj) = 
C^9l^^-l^^^Y{<l>+,cl>_,u), where r(</>+, </._,a.) ^ 

(tr(G+^(c.+ )<5^„G-Jc._)))^ and = 
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Apart from an overall factor w, this resembles the two- 
particle response function employed to compute the op- 
tical conduction properties of electronic systems. 

To make further progress, we describe the correlation 
function Y in terms of a low energy effective field the- 
ory. The technical details of this mapping [TT] are nearly 
identical to those of our earlier treatment of the one- 
dimensional rotor |4], and we here restrict ourselves to 
a brief sketch of the principal steps. We start from a 
representation of the correlation function Y in terms of 
Gaussian integral over superfields (8j: 



Here, P = E}^^ < 



E^l and P 



Here, the superfield ip — {4'N,x.a} where a = b, f 
distinguishes between commuting and anti-commuting 
components, and A = ± between retarded and ad- 
vanced components. The pre-exponential term is given 
by = V'A',+,bV'o,+,bV'o,-,bV'iv -,b and = 

diag |^G'^^(w+), G^^(a;_)^ is a matrix block-diagonal in 

advanced/retarded space. To make progress with this 
expression, we apply the color- flavor transformation |12j . 
an integral transform that trades the integral over ^ 
and Wo for the integration over an auxiliary field, Z: 
K{(j), uj)= J D{Z, Z) (. . . ) exp(-5[Z, Z]), where 

S[Z,Z] = -strln(l - ZZ) +strln(l - e'^^Ul zU^^Z), 

'str' is the supertrace ;8J, and we have temporarily sup- 
pressed the pre-exponential terms for notational simplic- 
ity. Here, Z = {Zj\ja,N'a'} is a bi-local supermatrix 
field, subject to the constraints Zb.b — Z^^, Zf f = 
— Zjf and |Zb.bZj!jj^| < 1. The anti-commuting blocks 

Za.a' and Zee,', a ^ a' are independent. Physically 
(cf. Ref. [H for a more extensive discussion), the field 
ZN,a;N',a' ~ 4'N,a,+4'N' ,a' ,- dcscribcs the pair propaga- 
tion of a retarded and an advanced single particle ampli- 
tude at a slight difference in frequency, ui, and Aharonov- 
Bohm flux tp = ^+ — (/)_. The structure of the action 
S[Z,Z] shows that at these values field configurations 
ZU^^ ^ Z near-stationary under the adjoint ac- 
tion of the Bloch-Floquet operator dominantly contribute 
to the field integral. The identification of these 'slow 
modes' is facilitated by passing to a Wigner representa- 
tion, Zni,N2 ^AT,*, where N = (A^i + N2)/2, and $ is 
dynamically conjugate to N. Due to the fast relaxation 
of the dynamics in the space of angular variables, the 
modes of lowest action Zjv_$ = Z^ depend only on iV. 
These angular zero modes then produce the low energy 
representation 

Y{ip,ij) = - J dN J dQe-^^^^ str(g^P)str(Qo^), 

S[Q] dNsti{D{d^Qf + 2iujQ<j\^). (3) 



ii^BP, where the 
2 X 2-matrices E^^^^-^-p act in the space of A/a-indices, 
carry a unity at position and zero elsewhere. The 

matrix field Q is given by 



Q = 



1 Z \ 3 / 1 Z 
Z 1 j ^AR, 1 z 1 



where cf\^ — -Bar ^ ^ar- The fiuctuations of these 
fields are governed by an action S[Q] identical to the 
action of the 'diffusive' nonlinear cr-model of disordered 
metals (subject to an Aharonov-Bohm fiux). The ro- 
tor's diffusion constant is given hy D — where 



C = (si'a.6dx\x=cos9VY is a constant of 0(1) whose 
detailed value depends on the kicking potential. (The 
perturbative integration over non-zero mode configu- 
rations Z^v^^ generates corrections to D of 0{\/K) 
which we do not discuss here.) Finally, 9^ = ((?„ -|- 
iip[a\^, ], 9„j , . . . 9„^_ J, is a covariant derivative ac- 
counting for the coupling to an AB flux in the compact 
ri-direction. 

Technically, Eq. ([3| represents the main result of the 
present paper. We have described the low energy physics 
of the quasi-periodic QKR in terms of the nonlinear a- 
modcl of disordered metals [13j. The construction is 
parametrically controlled by the parameters h/K,uj ^ 1 
and corrections to the effective action are small in these 
parameters. In the following, we discuss a number of 
physical predictions deriving from the representation (|3| . 

Metallic regimes. In dimensions d > 3 {q = oo) 
or d > 4 (g finite) the system supports an Ander- 
son (metal/insulator) transition. In the metallic phase, 
K/h 3> 1, fluctuations are weak and the action may be 
expanded to quadratic order in the generators Z. Doing 
the Gaussian integral over Z, one then obtains 



Y{ip,uj) 



1 



D^ip^ - ioj 



where ~ Z? is the diffusion constant weakly renor- 
malized by non-linear and frequency dependent correc- 
tions to the quadratic theory. Substituting this result 
into the expression for E{t), we obtain diffusive growth 
E{t) ^ Dt, corresponding to a finite optical conductivity. 

Localized regimes. For q = oo, the system is in a local- 
ized phase in low dimensions, d < 3, or below the Ander- 
son transition at K/h = 0(1) in d > 3. In these regimes, 
the diffusion constant is undergoing strong renormal- 
ization, '^^^ iuj. This in turn leads to saturation 
E{t) const., and vanishing (static) conductivity. For 
d > 3, the condition of scale invariance of the critical con- 
ductivity leads to the predictions {—iuj)^^ and 
E{t) ^ td at the critical point. This asymptotic agrees 
with the experimental observation [7] on the scaling of 
E{t) in d = 3. 
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Super-metallic regimes. For finite q < ^ and d < 3 
the system is localized in the virtual directions and de- 
localized along the real angular momentum direction (cf. 
Fig. [l] c)). Resonant transmission through the discrete 
levels of the ensuing system of effectively finite size then 
leads to 'super-metallic' growth E{t) = Ct^ at large time 
scales, and a corresponding diverging (static) conductiv- 
ity. The absence of rigorous descriptions of strong local- 
ization notwithstanding (for the exceptional case of d = 2 
see below) , we may apply phenomenological reasoning to 
estimate both the coefficient, C, and the crossover time, 
t^, to super- metallic scaling: at short times, the uncer- 
tainty in quasilevel resolution, ^ t^^, is larger than the 
characteristic quasilevel spacing Aj = l/{qLf^^) of a fic- 
titious system of size q x Lf^^, where Lt = {D^-ity^^ 
is the characteristic extension of a diffusive process of 
duration t in the virtual directions, and -Dt-i is a short- 
hand for the diffusion coefficient renormalized down to 
frequency scales w ^ t~^. The level mixing then leads 
to diffusive growth E{t) ^ Dt. The borderline condition 
= marks the crossover to long-time dynamics, 
t > t^, governed by localization effects. In this regime, 
individual levels are no longer mixed by quantum un- 
certainty. The coherent propagation through individual 
states then leads to E{t) = Ct^, where the coefficient 
C is fixed by the matching condition Dt^t^ = Ci^, i.e. 
C = DtJ^. 

The quantitative determination of both ^ and in- 
volves localization phenomena in slab geometries and is 
more difficult. Generally speaking, the above condition 
q < on the transverse extension refers to the bulk d- 
dimensional localization length. In contrast, the scale 
is determined by the localization length characterizing 
the quasi (d — l)-dimensional low energy regimes. The 
application of scaling arguments leads to ^ q^D 
and ~ exp(2q_D) in dimensions d = 2 and o? = 3, re- 
spectively. In dimensions d > 3, the situation is more 
complicated in that the system supports an Anderson 
transition in the underlying {d — l)-dimensional system. 
While the finite transverse extension of the system makes 
it difficult to determine the transition point, it is clear 
as a matter of principle that lowering the bare value of 
D ^ [K/Kf' will trigger a localization transition, which 
will manifest itself as a metal/super-metal transition in 
the long time scaling of the observable E. 

In view of the phenomenological nature of these ar- 
guments, it is reassuring that for d = 2 a much more 
sophisticated approach can be formulated. In this case, 
the theory is defined on an infinitely long cylinder, N — 
(n,ni) e [0,(7] X K. Fluctuations of the field Q inho- 
mogeneous in rt-direction are penalized by an action 
cost (cf. Eq. ^) of at least ^ D/q^, which we iden- 
tify as the inverse of the diffusion time in n-direction. 
For frequency scales smaller than that, n- fluctuations 
can be neglected, and we are left with the quasi one- 
dimensional field Qn = Qni with effective action S[Q] = 



-| / dni stT{D{dn,Q)'^+Dip^[Q,alT^]^+2iuj Qai^). The 
correlation function y ((^, w) of this theory can be com- 
puted by adaption [11] of Efetov's transfer matrix tech- 
nique [S] . The qualitative features mentioned above then 
follow from scaling properties of the corresponding solu- 
tions. For a quantitative discussion including pre-factors, 
we refer to Ref. [TT] . 

Summarizing, we have introduced a microscopic the- 
ory of quantum phase transitions in the quasi-periodic 
QKR. For irrational values of Planck's constant, the sys- 
tem is described by a d-dimensional nonlinear cr-model 
which entails a near perfect analogy to the physics of 
d-dimensional disordered metals. However, for rational 
values, its effective topology changes, and a dimensional 
reduction to a quasi {d — l)-dimensional system takes 
place. We discussed the ensuing consequences, includ- 
ing the existence of a metal/super-metal quantum phase 
transition in d > 4. It stands to reason, that the interplay 
of localization and these reduction phenomena should be 
observable in experiment. 
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